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1. Introduction 

Kaon physics is a rich resource when searching for new physics beyond the standard model. 
There have been a series of essential quantum jumps in the understanding of kaon spectrum, kaon 
decays, kaon mixing, partonic analysis of kaon, kaon form factors, and so on. In this article, we 
will review some of the recent progress in kaon physics. 

First, we review recent progress in understanding the pattern of taste symmetry breaking in 
the pion multiplet spectrum of staggered fermions. We compare a number of popular improvement 
schemes for staggered fermions and provide a guide to the best scheme by comparing the mass 
splittings of the pion multiplet in various improvement schemes. 

Second, we address the issue of calculating f n and f% with extremely high precision, which 
have a massive impact in determining some of the low energy constants such as L5 and the CKM 
matrix elements. We review the results of MILC (AsqTad staggered fermions in 2 + 1 flavor QCD), 
RBC/UKQCD (domain wall fermions in 2 + 1 flavor QCD), and NPLQCD (a mixed action scheme 
in 2 + 1 flavor QCD). 

Third, we review the phenomenological importance of indirect CP violation for the CKM 
unitarity triangle. We review the history of calculating Bk in quenched QCD, and compare the 
capabilities of the various improvement schemes for staggered fermions to reduce scaling violations 
in Bk- This comparison will guide us to the best scheme to improve staggered fermions. We 
will also review recent progress in calculating Bk in 2 + 1 flavor QCD using improved staggered 
fermions and domain wall fermions. 

Fourth, we address critical issues in calculating e'/e on the lattice. We explain a serious 
ambiguity in defining a lattice version of left-right QCD penguin operators, in particular Q(„ which 
was originally pointed out by Golterman and Pallante. We review a recent numerical study on 
this issue by RBC, which confirms the original prediction by Golterman and Peris. This leads 
to a conclusive claim that it is not possible to reliably calculate the leading contribution of to 
e'/e in quenched QCD. In addition, we also review theoretical progress in understanding the left- 
left QCD penguin operators by Golterman and Pallante. This explains the existence of another 
serious ambiguity in a lattice version of left-left QCD penguin operators, in particular Q\, Q2 and 
Q\, which leads to the conclusion that it is not possible to reliably calculate even the sub-leading 
contribution of Q4 to e'/e in quenched QCD. We also propose a solution to the problems in the 
end. 

Fifth, we review the history of calculating the hadronic matrix elements of K — > %% directly. 
We discuss the Maiani Testa no-go theorem, and consider a number of methods to get around it: the 
diagonalization method, the H parity boundary condition method, and the moving frame method. 
We also present recent numerical work in the moving frame by RBC. 

Sixth, we review the progress in calculating the vector form factor of decays on the lattice, 
which determines \V US \ with high accuracy. The results of RBC, HPQCD/FNAL, and JLQCD are 
presented. 

Seventh, we review the progress in calculating the first moment of the kaon distribution am- 
plitude. The results of UKQCD and QCDSF/UKQCD are presented. 

Finally, we summarize the whole discussion and end with some conclusions. In the final 
remarks, we propose to use Fat7 staggered fermions as the best candidate to generate 2 + 1 un- 
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quenched gauge configurations using the rational hybrid Monte Carlo algorithm (RHMC). 

2. Pion, Kaon spectrum 

Staggered fermions have four degenerate tastes by construction, which leads to SU(4) taste 
symmetry in the continuum limit a — ► 0. However, the taste symmetry is broken at a finite lattice 
spacing (a > 0). The splittings in the pion multiplet spectrum are a non-perturbative measure of 
the taste symmetry breaking. From a theoretical point of view, a natural tool to study the pion 
multiplet spectrum is staggered chiral perturbation theory |JT], Q). According to staggered chiral 
perturbation theory, the taste symmetry breaking in pion multiplet spectrum happens in two steps: 
at 0(a) « G(p ) the SU(4) taste symmetry is broken down to SO(4) taste symmetry, and at higher 
order of €?(a 2 p 2 ) the SO(4) taste symmetry is broken down to the discrete spin-taste symmetry SWa 
[jl|]. As a consequence, we expect that, at G(a 2 ) k, G(p 2 ), the pion spectra can be classified into 
5 irreducible representations of SO(4) taste symmetry: {75 (g)/}, {75 (g) £5}, {75 (g> ^}, {75 <8> ^5}, 
{75 (g) <^v}- The pion with £5 taste is the Goldstone mode corresponding to a conserved axial 
symmetry which is broken spontaneously. 

The splittings between pion multiplets reflect the taste symmetry breaking directly. As we 
improve the staggered fermion action and operators more, the splittings get smaller and the taste 
symmetry gets closer to SU(4). In other words, the splittings between pion multiplets can be used 
as a non-perturbative probe to measure how much of the taste symmetry is restored. 

There have been a number of proposals to improve staggered fermions: Fat7 AsqTad [Q], 
HYP [Q], Fat7 [Qj, and others. Unimproved staggered fermions have three serious problems: (1) 
large scaling violation, (2) large one-loop perturbative correction, and (3) large taste symmetry 
breaking, which are well established through a number of elaborate numerical studies. It turns out 
that HYP and Fat7 improvement schemes have the smallest one loop collections (<C 5 %) to bilinear 
operators [§], which resolves the problem of large perturbative corrections completely. The AsqTad 
and HYP improvement schemes are already being used extensively for numerical study of various 
physical observables, and it is important to determine which scheme is the best. Even though the 
perturbative comparison in Ref. [||] provides a hint that HYP and Fat7 schemes are better than the 
AsqTad scheme, it was not conclusive at that time due to a lack of non-perturbative evidence of 
improvement for the remaining two issues. In this context, it is quite important to compare the size 
of splittings between pion multiplets for various improvement schemes, which will tell us which 
improvement scheme is best to restore the SU(4) taste symmetry. We will address the issue of large 
scaling violation later when we talk about Bk ■ Here, let us focus on the issue of the taste symmetry 
breaking by looking into splittings among pion multiplets. 

In order to probe the size of the taste symmetry breaking for unimproved staggered fermions, 
we have calculated the pion multiplet spectrum in quenched QCD (j3 = 6.0 and 16 3 x 64 lattice) and 
preliminary results are presented in Ref. [[7|]. In Fig. [I], we show (am n ) as a function of quark mass 
am q for various pion multiplets. As discussed in Ref. [0], there are two independent methods to 
construct the bilinear operators using staggered fermions: (1) the Kluberg-Stern method and (2) the 
Golterman method. In Fig. [T[ the left (right) plot is obtained using the Kluberg-Stern (Golterman) 
method. Note that there is essentially no difference between the Kluberg-Stern and Golterman 
methods. We also observe that the splittings among pion multiplets are comparable to the light 
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Figure 1: (am n ) 2 vs. am q for unimproved staggered fermions: (left) using the Kluberg-Stern method; (right) 
using the Golterman method. 



pion masses, which implies that ff(a 2 ) « ff(p 2 ). In addition, note that the slopes are noticeably 
different for various tastes, which implies that &(a 2 p 2 ) terms are so large as to be noticeable in the 
case of unimproved staggered fermions, even though they are of higher order in staggered chiral 
perturbation theory. 

In Fig. ^, we compare the results of AsqTad staggered fermions with those of HYP staggered 
fermions. The AsqTad staggered fermion data comes originally from Ref. [JSj] whereas the HYP 
staggered fermion data comes from Ref. [Q]. Since AsqTad staggered fermions use a Symanzik 
improved gauge action with a = 0. 12fm and HYP staggered fermions use the unimproved Wilson 
plaquette action with a = O.lfm, we expect that the scaling violations associated with the back- 
ground gauge configurations are similar in both cases. In the case of AsqTad staggered fermions, 
the size of the taste symmetry breaking effect is comparable to the light pion masses, which implies 
that €?(a 2 ) « &(p 2 )- In addition, note that the slopes of pion multiplets with different tastes are the 
same, which implies that the 0'{a 2 p 2 ) terms are negligibly small for AsqTad staggered fermions. 

In the right-hand side of Fig. ||, we show data for HYP staggered fermions. In this plot, we 
observe two things. First, the splitting among pion multiplets are extremely suppressed within 2a 
of the data point error bar, which implies that Gia 2 ) -C 6{j) ). Second, we notice that the slopes of 
pion multiplets with different tastes coincide within statistical uncertainty, which implies that the 
ff(a 2 p 2 ) terms are so small that the effect is not detected at all. 

Therefore, from Fig. ||, we conclude that HYP staggered fermions are remarkably better than 
AsqTad staggered fermions with regards to reducing the taste symmetry breaking effect. From 
Fig. [I] and Fig. ^ the ff(a 2 p 2 ) terms are so large for unimproved staggered fermions that the effect 
was noticeable in the data, whereas these higher order terms are negligibly small in the case of 
AsqTad and HYP staggered fermions, which we can count as a big improvement. Therefore, we 
can rank the improvement efficiency of various staggered fermions with regards to taste symmetry 
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0.06 



Figure 2: {am n ) 2 vs. am q : (left) for AsqTad staggered fermions; (right) for HYP staggered fermions (using 
the Golterman method). 

restoration as follows: 

unimproved stag < AsqTad stag < HYP stag (2.1) 
More details on the comparison of various improvement schemes can be found in Ref. [Q]. 

3. f„ and f K 

A precise determination of low energy constants in QCD is one of the major goals that we 
want to achieve on the lattice. In particular, f n and fx are relatively cheap to calculate on the 
lattice and their analysis is straight forward. The real challenge is that we need to produce a large 
ensemble of gauge configurations which includes sea quark back-reaction (sea quark determinant) 
with Np = 2 + 1 flavors (m u = ^ m s ). This procedure is quite expensive computationally. Then, 
we compute valence quark propagators on the gluon backgrounds to calculate hadron spectrum 
and hadronic matrix elements. To get as much physics information as possible from the valuable 
gauge configurations with specific sea quark masses, we use many values of valence quark masses, 
i.e. not just the same values as the sea quark masses. This physical set-up is known as "partially 
quenched QCD". It has been studied with a number of lattice fermion formulations: (1) AsqTad 
staggered fermions, (2) Wilson clover fermions, (3) Wilson twisted mass fermions, (4) domain wall 
fermions, and (5) overlap fermions. In the following, we discuss the cases of AsqTad staggered 
fermions, domain wall fermions, and a mixed action approach (valence quarks are domain wall 
fermions and sea quarks are AsqTad staggered fermions). 

3.1 f n and fx using AsqTad staggered fermions 

Staggered fermions are fastest to simulate on the lattice, and have an exact axial symmetry 
which protects the quark mass from additive renormalization and makes it possible to calculate 
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Figure 3: f K vs. m x + m y (left) and m^/(m x + m y ) vs. m x + m y (right). The data are obtained using AsqTad 
staggered fermions in Nf = 2+1 partially quenched QCD. 



weak matrix elements reliably. The physical quark flavors are each represented by a staggered 
fermion field, so for three quark flavors there is an exact SU(3) flavor symmetry (when the quark 
masses are degenerate) as in the continuum, along with SU(4) taste symmetry of each of the stag- 
gered fermion fields (which is broken at finite lattice spacing but restored in the continuum limit). 
The hadronic operators can be taken to have a specific taste, so the lattice formulation with stag- 
gered fermions reproduces the continuum one at tree level in chiral perturbation theory. However, 
pions with all kinds of tastes appear in loops, and through this taste violations enter everywhere. 
Hence, it is necessary to take into account taste violations (directly related to discretization effects) 
within chiral perturbation theory in order to fit the data and extract physical information with high 
precision. A systematic tool to incorporate taste violations into the data analysis is the staggered 
chiral perturbation theory [III 0]. 

Another caveat regarding broken taste symmetry is that we take y/Det(D + m) as the fermion 
determinant to get a single taste per sea quark flavor. This is often called the "fourth-root trick". 
There is a possibility that y / Det(D + m) produces non-perturbatively violations of locality and 
universality in the continuum limit, which might cause staggered fermions not to reproduce QCD. 
However, recent arguments and numerical investigations give encouraging indications that this 
possibility does not occur. For more details on this issue, see Ref . [|T(]] . 

Recently, MILC has been generating dynamical gauge configurations with Nf =2 + 1 fla- 
vor sea quarks (AsqTad staggered fermions). These configurations include coarser lattices (a = 
0. 15fm), coarse lattices (a = 0. 12fm), fine lattices (a = 0.09fm) and super-fine lattices (a = 0.06fm). 
On these lattices, the strange sea quark mass (M s ) of simulation runs the range of 0.1m s < M s < 
1 .2m s . The light sea quark mass of simulation goes down to m u = 1 1 MeV (m n « 240 MeV). 
The physical volumes of the lattices range from (2.4/ra) 3 to (3.4/m) 3 . Using these gauge con- 
figurations, MILC has been calculating full pion multiplet spectrum and decay constants for the 
Goldstone pions and kaons. The results are presented in Fig. |3[ In this extensive calculation, MILC 
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observable 


value 


fx 
h 
hlfn 
V us 

m. v (2GeV,MS) 
m M (2GeV,MS) 
m s /m U)d 
L 5 


128.6 ± 0.4 ± 3.0 MeV 
155.3 ±0.4 ±3.1 MeV 
1.208(2)±? 4 
0.2223^ 
90(0)(5)(4)(0) MeV 
3.3(0)(2)(2)(0) MeV 
27.2(0)(4)(0)(0) 
2.0(3)(2) x 10" 3 



Table 1: Results of data analysis by MILC. 



concentrates on the analysis of masses and decay constants of Goldstone pions and kaons. This 
analysis also involves the masses of non-Goldstone pions and kaons; these enter in one-loop chiral 
logs at NLO in staggered chiral perturbation theory [|ll|]. In the data analysis, it is necessary to 
cut out the data points of heavier quark masses one by one in order to get good fits to the form 
suggested by staggered chiral perturbation theory. The most updated results are summarized in 
Table [jj. For more details, see Ref. [^]. 

3.2 f n and using domain wall fermions 

In the domain wall fermion formulation, the breaking of chiral symmetry is due only to the 
finite length of the 5th dimension, which allows a small additive renormalization to quark masses. 
We call this a "residual quark mass"; it is typically only a few MeV and hence under control. Do- 
main wall fermions are computationally typically two orders of magnitude more expensive than 
improved staggered fermions. Since we can use continuum chiral perturbation theory for the anal- 
ysis of domain wall fermion data, the theoretical interpretation is simpler and more straightforward 
compared to improved staggered fermions. 



parameter 


value 


lattice 


16 3 x 32 and 24 3 x 64 


U 


16 


aM 5 


1.8 


a 


0.125 fm 


am s 


0.04 


am Utd 


0.01,0.02, 0.03 



Table 2: Parameters of dynamical domain wall fermion simulation by RBC and UKQCD. 

The RBC and UKQCD collaborations use domain wall fermions with Iwasaki gluon action 
to generate dynamical gauge configurations with Np = 2 + 1 sea quark flavors. The simulation 
parameters are summarized in Table ||. These lattices are comparable with the coarse MILC lattices 
(a = 0.12 fm). These groups measured pion decay constants and pion spectra in partially quenched 
QCD, and the preliminary results are shown in Fig. |4[ They put elaborate effort into fitting the 
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Figure 4: af % vs. am va i (left) and (am n ) 2 / (m va i +m res ) vs. m va i (right). The data were obtained using 
domain wall fermions in Nf =2 + 1 partially quenched QCD. 



data of f n and m\/m q simultaneously to the NLO form suggested by partially quenched chiral 
perturbation theory. It turns out that the simultaneous fits fail to describe the data. They believe 
that this failure is due to the fact that the quark masses are too heavy to interpret the data based on 



chiral perturbation theory. For more details, see Ref. [|14j]. 

3.3 fx/ fn using a mixed action scheme 

Another possibility in partially quenched QCD is to use a mixed action scheme. This is the 



approach taken by the LHPC and NPLQCD collaborations []12p. In this approach, the sea quarks are 
AsqTad staggered fermions and the valence quarks are domain wall fermions. Specifically, they do 
HYP smearing over the MILC gauge configurations with Nf = 2 + 1 AsqTad staggered sea quarks 
and calculate the valence quark propagators using domain wall fermions over these configurations 
to measure fn/ fn- hi this approach, the valence quark masses are tuned such that the pion masses 
match the Goldstone pion masses of AsqTad sea quarks. 

In the continuum SU(3) chiral perturbation theory, one can show that fx / ' f% is given, at the 
next-to-leading (NLO) order, by 



h 
h 



1 +Xi og + %yL 5 (Li) (3-1) 



2 2 



y = K f2 71 (3-2) 

JK 

L 5 {f % ) = L 5 (ll) - Ij^^ifx/V) (3-3) 

where /I is a renormalization scale in chiral perturbation theory, %\og represents all the NLO chiral 
log corrections collectively, and fl is the physical value of pion decay constant. In Fig. |5[ we show 
Lsifn) calculated using the mixed action scheme. It is interesting to note that the final value of 
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Figure 5: Ls(f%) vs. m 2 n / f 2 calculated by NPLQCD, using a mixed action scheme. 
Ik I 'fx obtained by NPLQCD is consistent with MILC, even though they apply the continuum chiral 



perturbation theory to analyze data of the mixed action scheme. For more details, see Ref. Q12Q. 

4. Indirect CP violation and B K 

In nature, there are two kinds of CP violations in the neutral kaon system: indirect CP violation 
and direct CP violation. The neutral kaon states in nature are described as CP eigenstates with a 
tiny impurity which is parametrized as e. In experiment, e is determined precisely as follows: 

e = (2.280 ±0.013) x 10~ 3 xexp(/7i/4) (4.1) 

From the standard model, it is possible to express e in terms of K — K mixing as follows [ |I3| ] : 

e = C e exp(/7i/4) ImA, X B K (4.2) 
X = ReA c [T 7l 5 (x c ) - r] 3 S 3 (x c ,x t )] - Re^Tj 2 5 (^) (4.3) 
h = V*V id , Xi = mf/M^ (4.4) 

Q = (4.5) 

Here, note that 5; (the Inami-Nam functions) and T] ( are determined reliably from perturbation 
theory; they represent the remaining effect of heavy particles when they are integrated out. The 
kaon bag parameter Bk is defined as 

B = (Ko\W_ ~ 75 ~ Y5)d]\K Q ) 

K l(K \s Y ^ 5 d\0)(0\sY^5d\K Q ) 
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Figure 6: Bk determined from the unitarity triangle. 



B K = C{ll)B K {ll) (4.7) 
C(ji) = a s (n) m[l + a s ^)J 3 ] (4.8) 

where jl is a renormalization scale of the AS = 2 operator, and Bk is renormalization group invari- 
ant. Vij is the CKM matrix element of the standard model. The theoretical challenge is to calculate 
Bk with such a high precision that it over-constrains the unitarity of the CKM matrix. 

It is possible to determine Bk from the untarity triangle if we obtain the CKM matrix elements 
from other experimental results. In Fig. we show the result of Bk determined from the unitarity 
triangle. The result reported in CKM 2005 is B K = 0.68 ± 0.10. 



4.1 B K in quenched QCD 



There have been a number of efforts to calculate Bk in quenched QCD using staggered fermions 



Q15 , |lq ], twisted mass (TM) Wilson fermions [170, domain wall fermions [|18, |19p, and overlap 
fermions [20, 21]. In Fig. [7], the results of these calculations are summarized. We observe a large 
scaling violation in the calculation of Bk using unimproved staggered fermions by JLQCD [|l5|]. 
However, it is possible to reduce the scaling violation down to a negligibly small level using the 
HYP improvement scheme to improve the staggered fermion action and operators flT6[]. We will 



address the issue of scaling violation in detail in Sec. 4.2. An interesting demonstration of the 



progress in Bk calculations is that the value of Bk calculated using the TM Wilson fermions ends 



up being in agreement with other Bk calculations using different fermion formulations Q17|]. The 
calculation of Bk with TM Wilson fermions has been discussed in great detail in the plenary talk 
by Pena [22] so we skip all the details here. The values of Bk calculated by CP-PACS jl8| ] and 



by RBC [ ]19[ ] using domain wall fermions are in good agreement with those of other fermion for- 
mulations. The Bk values calculated by two independent groups [20, ^T] using overlap fermions 
are also in agreement with those of other fermion formulations, although they have relatively large 
error bars. In these calculations, it is assumed that m s = m^, i.e. the effects of non-degenerate quark 
masses are neglected. A simple-minded world average estimate of quenched Bk with degenerate 
quarks is fl^GeV.NDR) = 0.58(4). 
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Figure 7: B K {2GeV, NDR) calculated in quenched QCD. 



4.2 Scaling violation in Bk and its reduction with improved staggered fermions 

We now focus on the scaling violation in Bk calculated using staggered fermions. In Fig. ||, we 
show Bk data as a function of lattice spacing a (fm). Here, we observe that there is a large scaling 
violation (= discretization error) in the data calculated by JLQCD using the unimproved staggered 
fermion action and unimproved operators (constructed using thin links). The blue (green) data 
points in Fig. represent the Bk values calculated using the AsqTad (HYP) staggered fermion 



action and unimproved operators [J23p. We notice that in this case the scaling violations are reduced 
by about 50%. This comparison tells us that the HYP staggered fermion action is much more 
efficient in reducing the scaling violations than the AsqTad staggered fermion action. However, 
we also observe that the remaining scaling violations are still substantial. The red data points in 
Fig. ^represents the Bk values obtained using the HYP staggered fermion action and the improved 
operators with HYP fat links [jl^]. In this case, the Bk value at a = 0.1 fm is consistent with the 
world average which is obtained by extrapolation to the continuum limit of {a = 0). This implies 
that, once we improve the staggered fermion action and operators with HYP fat links, the scaling 
violations in Bk are reduced down to a negligibly small level. 

From Fig. || we conclude that we can rank the improvement efficiency of various staggered 
fermions with regards to reducing the scaling violations as follows: 

unimproved stag < AsqTad stag < HYP stag (4.9) 
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Figure 8: Z?R-(2GeV,NDR) calculated using various staggered fermion formulations. 

Combining this conclusion with that of Sec. ||, we see that HYP staggered fermions are remarkably 
better than AsqTad staggered fermions in reducing the taste symmetry breaking effect, the scaling 
violations, and the one-loop perturbative corrections simultaneously [Qj. 

4.3 B K in unquenched QCD with Np = 2 + 1 

Recently, there have been two serious attempts to calculate Bk in unquenched QCD with 
N F =2 + 1. In Fig. g, we show B K data in N F = 2 + 1 flavor QCD. On the left-hand side of 
Fig. |9], we show the Bk data of the RBC and UKQCD collaborations calculated using domain wall 



fermions at a = 0.125 fm [24]. On the right-hand side of Fig. BL we present the Bk data (a = 0.12 



fm) obtained using a mixed action scheme in which valence quarks are HYP staggered fermions 



and sea quarks are AsqTad staggered fermions ||25|]. The Bk data of domain wall fermions does 



not have a noticeable dependence on the non-degenerate quark masses. However, we observe that 
the Bk data of improved staggered fermions shows a noticeable dependence on the non-degenerate 
quark masses, corresponding to about 3 a difference between the degenerate and non-degenerate 
quark mass combinations near the physical kaon mass. Certainly, this is a puzzle. One possible 
explanation is that there is a significant difference in statistics between the two sets of data: 75 
measurements for the domain wall fermion data and 640 measurements for the improved staggered 
fermion data. In other words, it might be that this dependence on the non-degenerate masses only 
shows up with extremely high statistics. At any rate, this issue needs further investigation in the 
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aM q (avg) M K 2 (GeV 2 ) 

Figure 9: in unquenched QCD with Nf =2 + 1; (left) calculated by RBC and UKQCD collaborations 
using domain wall fermions and (right) calculated by Lee et al. using a mixed action scheme in which 
valence quarks are HYP staggered fermions and sea quarks are AsqTad staggered fermions. 



future. 

There has been another attempt to calculate Bk in unquenched QCD using the AsqTad stag- 
gered fermions and unimproved operators, where Bk is calculated with degenerate quark masses 
set to am s /2 p3|]. Since this work does not address the main issue of this section, we omit the 
details. 



5. Direct CP violation and Re(e'/e) 



Direct CP violation occurs in kaon decays via electroweak interactions where the CP sym- 
metry is explicitly broken. This process is parametrized by e'/e. Recently the KTeV and NA48 
collaborations have determined Re(e'/e) from high statistics experiments [^6], ^7|]. The results are 



Re(e'/e) = (20.7 ±2.8) x 1(T 4 (KTeV) 
= (14.7 ±2.2) x (NA48) 



(5.1) 
(5.2) 



Here, note that both KTeV and NA48 claim that Re(e'/e) is large and positive (»1.0x 10 4 ). 
From the standard model, we can derive the theoretical expression for e'/e as follows [13]: 



e'/e = lm(V*V td ) 



p(l/2)_ p (3/2) 



exp(z'<J>) 



10 



P {1/2) =r£ yi (M)<fi/(M)>o(l-QiB) 

i=3 
10 

P (3/2) = ^£ y . (M)((2 . (M))2 



(Q 



i 1 



r 

{nn(I)\Qi\K) 



(5.3) 
(5.4) 

(5.5) 
(5.6) 
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= G F (0 

r 2|£|Re(A ) 1 ' 

where n is a renormalization scale of the £2, operators. Here, f^ 1 / 2 ) and p( 3 / 2 ) represent the AI = 
1/2 and AI = 3/2 contributions to e'/e respectively. The Wilson coefficient functions y,-(ju) can be 
determined reliably from the perturbation theory in a specific renormalization scheme. Note that 
the phase shift <J> = <J> e / — <J> e and co = 0.045 (the AI =1/2 rule). The parameter £2fb represents 
the isospin breaking effect. 



As can be seen in Eq. 5.3, e'/e is obtained as a result of destructive interferences between 
pi 1 / 2 ) an d p( 3 / 2 ). In particular, p( 3 / 2 ) is fully dominated by electroweak penguin contributions, 
especially £2s- P' 1 / 2 ) on the other hand is governed by QCD penguin contributions, especially 
Q(,. The main sources of uncertainty in the calculation of e'/e are the hadronic matrix elements 
{nn{I)\Q,i\K) . Hence, the theoretical challenge on the lattice is to calculate {%Tt{I)\Qi\K) reliably. 

The RBC and CP-PACS collaborations reported that e'/e is small and negative. Their cal- 
culation is done in quenched QCD with operators defined in the SU(3) flavor group theory using 



domain wall fermions [|28|, |29|]. However, Golterman and Pallante have pointed out that there is an 



serious ambiguity in the lattice version of left-right QCD penguin operators (Qd and Qs) because 



the flavor group in quenched QCD is not SU(3) but SU(3I3) []30|]. It turns out that the ambiguity in 



<26 has such a large effect on e'/e that it can flip the sign of e'/e in quenched QCD: this has been 



demonstrated in calculations using HYP staggered fermions [31, B2J]. Recently, the RBC Collabo 



ration have investigated this issue in great detail [33], which will be discussed in the next section 



In addition, Golterman and Pallante point out that there is another ambiguity in Qi (i = 1,2,3,4) 
in quenched QCD, which has a serious effect on Re(Ao) and the sub-leading contribution of Q4 on 



e'/e [34]. In other words, it is not possible to calculate reliably even the sub-leading contribution 
of Q4 as well as the leading contribution of to e'/e in quenched QCD. We will address this issue 
later. 



5.1 Ambiguity in £?6 and Q5 

We now focus on the left-right QCD penguin operators, especially which makes the domi- 
nant contribution to pt 1 / 2 ). In the standard model, Q(, is defined as 

Q6 = sr^(i-r 5 )dY,qM^ + r5)q (5.8) 
i 

£26 belong to the (8,1) irreducible representation of the SU (3)l SU (3)r flavor symmetry group. 
Unfortunately, £26 mixes with £2,- (z = 3,4,5,6) in (partially) quenched QCD, where for example, 
£26 is defined as 

26 = ^(1" 75^X^(1 + 75>7 (5-9) 
q 

where q represents a ghost quark. Here, note that the flavor symmetry of (partially) quenched QCD 
is not SU (3) but SU (3 + Np\3) where Np is the number of sea quark flavors. In quenched QCD 
(Nf = 0), Golterman and Pallante propose to decompose £26 as follows: 

Qe = Qf + Qf s (5-10) 
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Figure 10: in quenched QCD calculated by the RBC collaboration. 



Qf = ^(\-y 5 )dY J Qy^i + 75)Q 



a 



QNS 



(5.11) 

(5.12) 
(5.13) 



where Q = (u,d,s,u,d,s) and A = diag(l, 1,1,-1,-1,-1) is a diagonal matrix. See Ref. [30, 33] 
for the details of the group theoretical interpretation on this decomposition. s can be represented 
by a set of operators in quenched chiral perturbation theory which have analogues in continuum 
chiral perturbation theory. However, Q® N introduces a new low energy constant which is called 



CCq*. In addition, Golterman and Peris have predicted that « 600^'^ J35|], which implies that 
the matrix elements of Q(, are completely dominated by unphysical leading contribution from the 



lattice artifact off s . Unfortunately, the work by the RBC and CP-PACS collaborations in Ref. |^8 



2% contains this artifact as it is. As mentioned, the work of Ref. [|31|, |32|] shows clearly that this 
ambiguity is so large that it can even flip the sign of e'/e. 

Recently, the RBC collaboration has measured the low energy constant a^ s by two indepen- 
dent methods [|33|]: (1) using the hadron matrix element of K — > which is originally suggested 
by Golterman and Pallante, and (2) using the K + — > n + matrix element. In the case of domain 
wall fermions, the 2nd method has the advantage that one does not have to worry about the residual 
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quark mass m res , because the divergent NLO term vanishes in the limit of degenerate quark masses. 
In Fig. [k], we show the procedure to obtain a^ s by extrapolation to the chiral limit. The figure 
demonstrates that both methods produce the same result in the chiral limit and it is consistent with 



the prediction by Golterman and Peris. See Ref. [ ]33[ ] for more details. 

5.2 Ambiguity in <23 and Q4 

Here, we focus on the left-left QCD penguin operators (<23, Q4, part of Q\ and Q2, • • •)■ As an 
example, let us consider Q3 and Q4: 

Q3 = Sa7[i{l ~ 75)d a Y,qp7[i(l ~ 7s)qp = (s a da) L ^P ( }p)L (5-14) 
<z q 

Q4 = s a y^(l ~ 75)dp^qp7Li( l ~ 7s)q a = (s a dp) L Y,(vp1cc)l (5-15) 

q q 

Both £?3 and Q4 belong to the (8,1) irreducible representation of the SU (3)l <£> SU (3)r flavor sym- 
metry group, which also corresponds to / = 1/2. For the following discussion it is convenient to 
change the operator basis as follows: 

Q± = Qs±Q4 (5.16) 

Here, note that Q (<2+) is symmetric (anti-symmetric) in the two anti-quark flavor indices and 
two quark flavor indices. Once more, the problem is that the flavor symmetry group of (partially) 
quenched QCD is not St/ (3) but St/ (3 +N F \3). Hence, Q± is reducible under St/ (3 +N F \3). We 
can decompose Q± as follows: 

e± = ^ei+ei (5.n) 

Q± = {qai^qa)L{qpqp)L±{q a ^qp)L{qpqa)L (5.18) 

21 = (q a Aq a ) L (qpAqp l ) L ±(q a Aqp) L (qpAq a ) L (5.19) 

where q represents (valence plus sea) quark fields as well as ghost quark fields. The spurion fields 
A and A are defined as 

Ay = 8 i3 8 j2 (5.20) 

3 3 3 3 3 

A = diag(l ,1 ,1 , ,) (5.21) 

BK N F N f N f N f N f ' 

In this decomposition, the Q± part produces those terms which have analogues in continuum SU(3) 
chiral perturbation theory. However, the Q± part introduces new low energy constants which are 



pure lattice artifacts. In quenched QCD, the decomposition is slightly different; see Ref. [ |34| ] for 
the details. At any rate, this unphysical lattice artifact from the <2± part appears at the leading order 
in (partially) quenched chiral perturbation theory. This has a serious effect on Re(Ao) through Q\ 
and Q2 as well as e' / e through Q4. Consequently, it is not possible to calculate reliably not only the 
leading contribution of (?6 but also the sub-leading contribution of Q4 to e'/e in quenched QCD. 



Moreover, the analysis of Re(Ao) given in Ref. [28, 29, Bq] did not take into account the effect of 



this lattice artifact from the <2± part, which implies that the results are contaminated by this as well. 
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5.3 Solution to the problem 



In the previous sections |5T| and p.2| , we explained that it is not possible to reliably calculate 
e' /e in quenched QCD due to a fundamental ambiguity in and Q4. Is there any way to get 
around the problem? One possible solution is the following 

• Do the numerical study in partially quenched QCD with Np = 3 or Nf = 2 + 1. This insures 
that the low energy constants are physical. 

• Calculate the hadronic matrix elements of and Q5 using an singlet operator of SU (6\3)r 
to remove the unphysical lattice artifact. 

• Calculate the hadronic matrix elements of Q4 and (?3 using linear combinations of the Q± 



operators defined in Eq. 5.18. 



In this way, the calculations will extract only those low energy constants which are physically 
relevant to the direct CP violation without contamination by non-physical lattice artifacts. 

6. K — > tztz decay 

There have been two independent approaches to calculating the non-leptonic kaon decay pro- 
cess K — > 7171 on the lattice: a direct approach and an indirect one. In the indirect approach, we 
calculate the hadronic matrix elements of K — > % and K — > and reconstruct K —> Tin amplitudes 
out of them using the chiral perturbation theory. Since it is relatively easy and computationally 
cheap, the indirect method has been widely used to calculate the K — > Tin amplitudes [28, 29], pi] ]. 



In the direct approach, we calculate the K ^ tctc matrix elements directly with the final state pi- 
ons carrying a physical momentum of \p\ = 206 MeV. This method has a number of difficulties 
on the lattice. In Ref. [[37|], Maiani and Testa pointed out a "no-go theorem": we can not obtain 
K — ► 7t(p)7l(—p) but only K — ► 7r(0)7r(— 0) on the lattice. This is seen as follows. We want to 
obtain the matrix element of (n(p)n(—p)\Q\K) by calculating the following four-point function 
G: 

G{t,t n ,t K ) = (o\Mp,t*)M-p,tMt)MtK)\o) (6.1) 

However, the asymptotic behavior of G is 

lim G(f ) W J 0=Z- r exp(-£ _ g (^-f)) • K0)^(0)|Q|Z) -Z^-exp^^-^)) (6.2) 
where 

Z K = (K\(j) K \0) (6.3) 

z# = (o\Mp)M-p)\x(qM-q)) (6.4) 



Emcjt = 2 \fm 2 Tc+P 2 (6-5) 

Hence, we need {%{p)%{— p)\Q\K), but we can get only the information on (7r(0)7r(0)|<2|^0, be- 
cause 1 7r(0)7r(0)) is the ground state of two pions with periodic boundary condition in the spatial 
direction. 
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6.1 Circumventing the Maiani-Testa theorem 

Is there any way to get around the Maiani-Testa theorem? Luscher and Wolff proposed the di- 



agonalization method in Ref. [38]. In this method, we calculate the following correlation function: 



c nm {t) = (n„(*)n m (o)) (6.6) 

&n{t) = $x{j>n,t)ty%{-Pn,t) (6.7) 

Then, we can construct the following matrix M(t,to): 

M(t,t ) = C(? r 1/2 • C(t) ■ C(f r 1/2 (6.8) 

By diagonalizing the M(t,to) matrix, we can obtain the two pion eigenstates and eigenvalues. Of 
course, this takes enormous amount of computing resources because we need to take a large number 
of momenta in order to get a few low eigenvalues, and the signal gets poor as we insert the pion 
momentum. A number of lattice groups applied this method to the 1 = 2 elastic pion scattering 



problem [39, 4 



One way to get around the Maiani-Testa theorem is to modify the boundary condition (BC) 



for the pions [41]. Of course, there are many different choices to modify BC: (1) H parity BC, (2) 



G parity BC, (3) twisted BC and so on. As an example, let us consider the H parity BC. 

H(u,d) = (-u,d) (6.9) 
H|w ± ) = -|w ± ) (6.10) 

H\7l°) =+|7T°) (6.11) 

The H parity BC causes the ground state of n ± to have a finite momentum \p\ = \/2(2n/L). 
Hence, if we construct two pion state out of n + and %~ with the H parity BC, we can calculate 

(n(p)n(-p)\Q\K). 

Another way to get around the Maiani-Testa theorem is to use the moving frame (LAB) in- 



stead of the center of mass (CM) frame. In Ref. [|42|, [43| [44]], it was proposed to put a non-zero 
total momentum into the final two pion state and initial kaon state. In other words, calculate 
(n(P)7t(0)\Q\K(P)) . Of course, we need a formula which can convert this result on the LAB 
frame into that on the CM frame. 

6.2 Lellouch-Luscher formula in the moving frame 

Recently, Lellouch and Luscher have provided a formula in the CM frame which connects the 
K ^ Tin amplitude calculated in a finite box (on the lattice) with the K — ► %% amplitude in the 



infinite volume [45]. Let us define A as the on-shell decay amplitude in the infinite volume of the 



CM frame and M as that in the finite volume of the CM frame. Their relation is 

| A p = 8l (^) 3 { p £|M + ,^)}| M | 2 (6 . 12 ) 

where 8 is the scattering phase shift and is obtained from the following relationship: 



2\ ml + p 2 = m K (6.13) 
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q = Yn P 
d{p)=rni-<^(q) 



tan(<%)) 

Zoo(i;<? 2 ) 



3/2 



Zoo(i;# 2 ) 
i ^ i 



(6.14) 
(6.15) 

(6.16) 
(6.17) 



However, this formula is valid only in the CM frame. Hence, we need a similar relationship which 
works in the LAB frame. The latter is given in Ref. p2| , [43| , Q. We can use the Lellouch-Luscher 
formula with the following modifications for K(P) — > 7r(/ , )7r(0): 



(^ b ) 2 =^L+p 2 



7 



~|M C m| 



pLAB 



(6.18) 
(6.19) 

(6.20) 



The scattering phase shift is modified as follows: 



2 \ m K + P = m K 



q= Yn P 

d{p)=nn-§p{q) 



tan (^(4)) 

Zoo(i;^ 2 ) = 



JqTt 



3/2 



4?r "i + «2 + r 2 («3 + LP/47T) 2 - <? 2 



where we assume that P = (0,0,27r/L). The relation between |A| and |M£ab| is 

5 < d8(p) , d$p(q) 



i ,i2 o 2-( E nn\ f d8(p) c 
\A\ =**?[—) {P^ + «- 



dq 



M. 



LAB\ 



(6.21) 

(6.22) 
(6.23) 

(6.24) 
(6.25) 



(6.26) 



Note that Eq. 6.26 becomes identical to Eq. 5.12 in the limit of P = and 7=1 



6.3 Numerical study in the moving frame 

Recently, the RBC collaboration has calculated the K(P) — ► 7r(P)7r(0) amplitudes in the LAB 
frame in quenched QCD [46]. This numerical study is performed at \ ja = 1.3 GeV on the 16 3 x 32 
lattice with L s = 12 using domain wall fermions and the DBW2 gauge action. Using the formula of 
Eq. |6.26 , the result in the LAB frame is converted into that of the infinite volume in the CM frame. 
A preliminary result for Re(A2) is shown in Fig. O. 
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Figure 11: Re(A2) vs. m^. in quenched QCD calculated using domain wall fermions. 

7. K i3 decay 

The £/3 decays are the kaon beta decay processes: K + — > 7T°Z + V/ and /T — > 7T~/ + V/, where Z 
represents the electron or the muon. These decay channels play an important role in determining 
the CKM matrix element \V IIS \. The decay rate of the processes is given in Ref. [07j]: 



r = ^ i4-|v»| 2 -c 2 • |/ + (0)| 2 •/ • (i + 8) 



(7.1) 



where C is a Clebsch-Gordon coefficient (1 for K° decay; l/y/2 for K + decay), / represents the 
phase space integral, and 8 represents radiative corrections from electroweak and electromagnetic 
interactions. Here, /+(0) is a form factor which is defined for the neutral kaon decay as 

{%{p')\sy^u\K(p)) = (p„ +p' tl )f + {q 2 )+q^f-{q 2 ) (7-2) 

where q = p — p' represents the momentum transfer. For later discussion, it is convenient to define 
the scalar form factor fo{q 2 ) as 



Hq 2 )=U(q 2 ) + ^— 2 



f-(q 2 



_ /-(<? 2 
" U(q 2 



(7.3) 
(7-4) 
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Figure 12: fo(q 2 ) in two flavor QCD calculated by the RBC collaboration. 



The decay rate T is well determined from experiment. From Ref. [48], 



1^/+ (0)| =0.2173 (7.5) 

Hence, if we know /+ (0) with extremely high accuracy, we can determine V us with as much pre- 
cision; this can impose a stringent constraint on the unitarity of the CKM matrix. The theoretical 
challenge on the lattice is to determine /+(0) with such a high precision that the uncertainty is 
under control within <1%, 

Vector current conservation implies that /+(0) = 1 in the exact SU(3) limit of m u = ma = m s . 
The Ademello-Gatto theorem implies that the deviation of /+(0) from unity is quadratic in the 
quark mass difference (m s — m u ). If we expand /+(0) in chiral perturbation theory [p9|], 

/+(0) = l+/2 + /4 + ^(p 6 ) (7-6) 

where /; represents terms of order 0(p l ). Here, f% does not depend on any new low energy con- 
stants which can be determined in terms of pseudoscalar meson masses and decay constants. How- 
ever, /4 does depends on new low energy constants. 



The lattice technique to calculate /+(0) on the lattice is the double ratio method jpOQ. 

{Ti\sy A u\K){K\uy A s\n) (m K + m x f 2 2 

{n\uy A u\%){K\sy A s\K) 4m K m n lJo{%>1 1 " ; 

where the pion and kaon states carry zero momenta and = (mx — m n ) 2 . Here, note that /+(0) = 
/o(0). 

Eq. ^7| determines /o(#o)> but we need to know /o(0) = /+(0). Hence, it is necessary to 
either extrapolate or interpolate in q 2 to the (Minkowski) q 2 = point at fixed pseudoscalar meson 
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collaboration 


/+(0) 


RBC 


0.968 (9) (6) (N F 


= 2) 


HPQCD/FNAL 


0.962(6) (9) (N F = 


2 + 1) 


JLQCD 


0.952(6)(-) (N F 


= 2) 


Leutwyler 


0.961(8) 





Table 3: Results for /+(0). 



masses. In other words, we need to insert small momenta in the kaon and pion states to change 
q 2 = (p — p') 2 (defined in Minkowski space) from q\ = (nix — m n ) 2 to negative q 2 such that we 
can interpolate fo(q 2 ) to the q 2 = point. For this purpose, we can calculate the following ratio on 
the lattice: 

m K + m n {n(p')\sy A u\K(p)) 



F(p,p') 



f Q (q 2 )\ E K (j>)+E K (j>i) 

where p and p' are the kaon and pion momenta respectively, Ek{p) is the energy of kaon with a 
momentum p, and E n {p') is the energy of pion with a momentum p'. Here, note that F(0,0) = 1. 



E K {p)+E n {p') (n(d)\sy 4 u\K(d)) 

Mf)f 1 + E K W-E x tf) ^ (78) 



In order to extract fo(q ) from Eq. we first need to know £ (q ). It can be directly determined 
by measuring the following ratio: 



R(pJ) 



{n{p')\sYAK{p)){K{p')\sy A s\K{p)) 

(7.9) 



{K{p')\sy A u\K{p)){K{p ! )\sy iS \K{p)) 

[(p + P')i + (p - pU (q 2 )] ■ [Ek(p) + E K (p')} 



l(E K (p) +E n {p')) + (E K (p) -E n {p'mq 2 )] ■ [(p + p'U 
After determining %(q 2 ) from this, it is then possible to determine fo(q 2 ) by combining Eq. 



7.8 and 7.9 



Recently, a main goal on the lattice has been to calculate /o(0) in unquenched QCD. Three 
groups have reported results on this. First, the RBC collaboration have calculated fo(q 2 ) in two 
flavor QCD using domain wall fermions at a = 0.12 fm on the 16 x 32 lattice [Mh. In Fig. O, we 



show their results for fo(q 2 ) as a function of (Minkowski) q 2 . Second, the JLQCD collaboration 
has calculated fo(q 2 ) in two flavor QCD using non-perturbatively G{a) improved Wilson fermions 



[p2[]. In this calculation, the Wilson fermion action is G\d) improved but the vector current operator 
is not fully G(a) improved. Third, the HPQCD/FNAL collaboration has calculated fo{q 2 ) in 2 + 1 



flavor QCD [p3[]. They use unimproved staggered fermions for the u and d quarks and an improved 
Wilson fermion for the s quark. They calculate only /o(#o) using the double ratio in Eq. 1JT_ and 



extrapolate to q 2 = using the pole model by setting the free coefficient to the experimental value. 
These results are summarized in Table |3[ 

8. Kaon distribution amplitude 

Exclusive reactions with specific hadrons in the final and initial states have received a lot of 
attention. The reason is that they are dominated by rare configurations of the hadron's constituents. 



22 



Progress in Kaon Physics on the Lattice 



Weonjong Lee 



One of these rare configurations is based on the hard mechanism in which only valence quark 
configurations contribute and all quarks have small transverse separation. The structure of the hard 
mechanism is relatively simple to understand. These hard contributions can be calculated in terms 
of hadron distribution amplitudes (DA's) which describe the momentum fraction distribution of 
partons at zero transverse separation with a fixed number of constituents. In particular, the leading- 
twist DA's of the pion and the nucleon have been the focus of much attention. The DA's of the 
kaon are important for understanding the exclusive B decays such as B — ► KK* in the framework 
of QCD factorization, perturbative QCD and light cone sum rules. 

The theoretical description of DA's for the kaon is based on the relation: 

(0\u(z)r v r 5 P(z,-zX-z)\K + (p)) z 2 =0 = ifKPvf d£ exp(i^-z) fa{^^ 2 ) (8-1) 

where z. v is a light-like vector, z 2 = 0, jl is a renormalization scale, and P(z, —z) is a path-ordered 
Wilson line connecting the quark and the anti-quark fields in a gauge invariant way: 

P(z, -z) = & exp ( - ig jP dy^ (y)j (8.2) 

Here, fa is the leading twist DA of the kaon. The physical interpretation of E, is that the s quark 
carries an .x = (1 + £)/2 fraction of the kaon momentum and the u anti-quark carries the rest (a 
1— jc = (1 — £)/2 fraction). The fa parametrizes the overlap of the kaon state with longitudinal 
momentum p with the lowest Fock state consisting of a quark and an anti-quark carrying the mo- 
mentum fraction xp and (1 — x)p, respectively. In Eq. |8.1| , the convention for the normalization of 
fa is as follows: 



1 dt<t> K (Z,n 2 ) = l (8.3) 
-l 



The conformal expansion provides a convenient tool to separate the transverse modes and the 

longitudinal modes. The transverse mode is formulated as a scale dependence of relevant operators 

3 /2 

and the longitudinal mode is described in terms of Gegenbauer polynomials C„ (£,)'. 

fa($,LL 2 ) = 1(1- (l + t «n(M 2 )d /2 (S)) (8.4) 



We refer to the a„'s as "Gegenbauer moments". 

Applying the light-cone operator product expansion (OPE) to the left-hand side of Eq. 8. 1 
can drive the following relationship for the n-th moment (<!;") of the kaon DA: 



we 



<nGu 2 ) = f^^fa^f) (8.5) 
, ,u. = i"u(0)Y, Y5 D m - D, n s(0) (8.6) 

where the operator is defined in Minkowski space, D=D — D, and {•••} denotes the traceless 
symmetrization of all the indices. The moments are related to the Gegenbauer moments a n as 
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Figure 13: (%)k in 2 + 1 flavor QCD calculated by the UKQCD collaboration using domain wall fermions. 



follows: 

ax = 5 -{£,) (8.8) 
a 2 = l(5(^ 2 )-\) (8.9) 

The theoretical mission on the lattice is to calculate the DA moments (^") in Euclidean space. 
The first step is to transcribe an operator in Euclidean space as follows: 

= B M1 • • ■ D M „ *(0) (8.10) 

Since the symmetry group on the lattice is the finite group H(4), which is a subgroup of 0(4), we 
need to classify On^ ^ . » in terms of irreducible representations of H(4). The latter, in general, 
contains a number of irreducible representations of 0(4). For the first moment of the kaon, there 
are two types of operators: 



r 4i = + 



(8.11) 



The operator requires a nonzero momentum component in the 1 -direction. The 0^ operator 
can be evaluated at zero momentum. 

Recently, the UKQCD collaboration has calculated the first moment of the kaon DA in 2 + 1 



flavor QCD using domain wall fermions [p4|]. They use the operator of type. In Fig. 13, we 
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Figure 14: (%)g in two flavor QCD calculated by the QCDSF/UKQCD collaboration using ff(a) improved 
Wilson fermions. Top: using Bottom: using ff^. 



show (%)k as a function of m s — m u d. Chiral perturbation theory predicts that there is no correction 
with chiral logarithms [55]. The data is certainly consistent with this prediction. 

The QCDSF/UKQCD collaborations have also calculated the first moment of kaon DA in two 
flavor QCD using &{a) improved Wilson fermions [56]. In Fig. 14, we show (%)k as a function of 
m K ~ n ^k- m this calculation, they use operators of both type and type. 



9. Summary and conclusion 

In this paper, we have reviewed recent progress in calculating kaon spectrum, Bk, s' /e, K — > 
7171 amplitudes, kaon semileptonic form factors, and moments of the kaon distribution amplitude 
on the lattice. Especially, there has been much progress in calculating physical observables in 
unquenched QCD using various lattice fermion formulations. For example, in the case of Bk, the 
quenched results of various groups now agree with each other within statistical and systematical 
uncertainty. 

Through elaborate numerical studies on the pseudoscalar meson multiplet spectrum and Bk, 
we learned that HYP staggered fermions are remarkably better than AsqTad staggered fermions 
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in reducing both taste symmetry breaking effects and scaling violations. However, HYP fat links 
are quite expensive when we calculate the fermionic force term in the rational hybrid Monte Carlo 
(RHMC) algorithm, because there are three SU(3) projection steps for each HYP fat link. For this 
reason, in Ref. []|], we proposed the Fat7 scheme to improve staggered fermions for the unquenched 
QCD simulation. The Fat7 scheme shares the same advantage as the HYP scheme in that the one- 
loop corrections for general composite operators are identical in both schemes In addition, the 
Fat7 scheme is computationally an order of magnitude cheaper for the the unquenched QCD sim- 
ulation using the RHMC algorithm, mainly because we may have only one SU(3) projection step 
for each Fat7 link and the memory usage is only a quarter of that for the HYP scheme. Taking into 
account all these facts, we recommend using the Fat7 staggered fermions to generate unquenched 
gauge configurations using the RHMC algorithm in the future. 

There are a number of fundamental difficulties for finding a good lattice version of the and 
£?4 operators in quenched QCD. A solution to this particular problem has been proposed here. 

In the case of direct calculation of K — > %% amplitudes, there have been a number of diffi- 
culties, in particular the Maiani Testa no-go theorem. A number of methods to get around the 
latter have been proposed and tested extensively. At present, the moving frame method looks quite 
promising and deserves further investigation. The twisted boundary condition method also looks 
promising but has not been discussed extensively here due to lack of space. It also deserves further 
attention in the future. 

Most of the topics in kaon physics have their chiral behavior known from (partially) quenched 
chiral perturbation theory. However, the corresponding results in staggered chiral perturbation 
theory are not available in many cases. This also requires our attention and efforts in the future. 

10. Acknowledgments 

Proof-reading by David Adams is acknowledged with gratitude. Helpful discussion with 
S. Sharpe is also acknowledged with many thanks. W. Lee acknowledges with gratitude that the 
research at Seoul National University is supported by the KOSEF grant (R01-2003-000-10229-0), 
by the KOSEF grant of international cooperative research program, by the BK21 program, and by 
the DOE SciDAC-2 program. 

References 

[1] Weonjong Lee and Stephen Sharpe, Partial Flavor Symmetry Restoration for Chiral Staggered 
Fermions, Phys. Rev. D60(1999) 094503, [hep-lat/9905023]. 

[2] C. Aubin and C. Bernard, Pion and kaon masses in staggered chiral perturbation theory, Phys. Rev. 
D68 (2003) 034014, [hep - 1 at / 3 4 1 4 ] . 

[3] K. Orginos, D. Toussaint, and R.L. Sugar, Variants of fattening and flavor symmetry restoration, 
Phys. Rev. D60 (1999) 054503, [hep - 1 at / 9 9 3 3 2 ] . 

[4] A. Hasenfratz and F. Knechtli, Flavor symmetry and static potential with hypercubic blocking, 
Phys. Rev. D64 (200 1 ) 034504, [hep-lat/0103029]. 

[5] Weonjong Lee, Perturbative improvement of staggered fermion using fat links, Phys. Rev. D66 (2002) 
114504, [hep-lat/0208032]. 



26 



Progress in Kaon Physics on the Lattice 



Weonjong Lee 



[6] Weonjong Lee and Stephen Sharpe, One-loop Matching coefficients for improved staggered bilinears, 
Phys. Rev. D66 (2002) 114501, [hep-lat/0208018]. 

[7] T. Bae, et al, Pion spectrum using improved staggered fermions, PoS LAT2006 (2006) 166, 

[hep-lat/0610056]. 

[8] C. Aubin, et al., Light pseudoscalar decay constants, quark masses, and low energy constants from 
three flavor lattice QCD, Phys. Rev. D70 (2004) 114501, [hep-lat/0407028]. 

[9] MILC collaboration, C. Bernard et al, Update on the physics of light pseudoscalar mesons, PoS 
LAT2006 (2006) 163, [hep-lat/0609053]. 

[10] S. Sharpe, Rooted staggered fermions: good, bad, or ugly?, PoS LAT2006 (2006) 022. 

[11] C. Aubin and C. Bernard, Pseudoscalar decay constants in staggered chiral perturbation theory, 
Phys. Rev. D68 (2003) 0740 11, [hep-lat/0306026]. 

[12] S.R. Beane et al., fx./ fn in full QCD with domain wall valence quarks, [hep- 1 at / 60 602 3]. 

[13] A.J. Buras, Flavor Dynamics: CP Violation and Rare Decays, [hep-ph/ 010133 6]. 

[14] RBC and UKQCD Collaborations, Meifeng Lin et al, Light pseudoscalar decay constants from 2 + 1 
flavor domain wallfermion simulation, PoS LAT2006 (2006) 185. 

[15] JLQCD Collaboration, S. Aoki et al., Kaon b parameter from quenched lattice QCD, Phys. Rev. Lett. 
80 ( 1 998) 527 1 -5274, [hep- 1 at / 9 7 1 7 3] . 

[16] W. Lee et al., Testing improved staggered fermions with m s and Bk, Phys. Rev. D71 (2005) 094501, 

[hep-lat/0409047]. 

[17] Alpha Collaboration, P. Dimopoulos et al., A precise determination ofB^ in quenched QCD, 
Nucl. Phys. B749 (2006) 69-108, [hep-lat/0601002]. 

[18] CP-PACS Collaboration, A. Ali Khan et al, Kaon B parameter from quenched domain wall QCD, 
Phys. Rev. D64 (2001) 114506, [hep-lat/0 105020]. 

[19] RBC Collaboration, Y. Aoki et al, The Kaon B-parameter from quenched domain-wall QCD, 
Phys. Rev. D64 (2001) 114506, [hep-lat/0 105020]. 

[20] MILC Collaboration, T. DeGrand, Kaon B parameter in quenched QCD, Phys. Rev. D69 (2004) 
014504, [hep-lat/0309026]. 

[21] R. Babich et al, K° — K° mixing beyond the standard model and CP-violating electroweak penguins 
in quenched QCD with exact chiral symmetry, [hep- 1 at / 605016]. 

[22] C. Pena, Twisted mass QCD for weak matrix elements, PoS LAT2006 (2006) 019. 

[23] HPQCD/UKQCD collaborations, E. Gamiz et al, Unquenched determination of the kaon parameter 
B^from improved staggered fermions, Phys. Rev. D73 (2006) 114502, [hep-lat/0603023]. 

[24] RBC and UKQCD collaborations, S.D. Cohen and J. Noaki, Bk on 2 + 1 flavor Iwasaki DWF lattices, 
PoS LAT2006 (2006) 080. 

[25] J. Kim, et al, Bk in unquenched QCD using improved staggered fermions, PoS LAT2006 (2006) 086, 
[hep-lat/06 10057]. 

[26] KTeV Collaboration, Measurements of Direct CP Violation, CPT Symmetry, and Other Parameters in 
the Neutral Kaon System, Phys. Rev. D67 (2003) 012005, [hep-ex/0208007]. 



27 



Progress in Kaon Physics on the Lattice 



Weonjong Lee 



[27] NA48 Collaboration, Direct CP violation in neutral kaon decays, Pramana 62 (2004) 601-606, 
[hep-ex/0303037]. 

[28] RBC Collaboration, Kaon matrix elements and CP violation from quenched lattice QCD: 1. The 
three flavor case, Phys. Rev. D68 (2003) 114506, [hep-lat/01 10075]. 

[29] CP-PACS Collaboration, Calculation of non-leptonic kaon decay amplitudes from K — > % matrix 
elements in quenched domain-wall QCD, Phys. Rev. D68 (2003) 114501, [hep-lat/0108013]. 

[30] M. Golterman and E. Pallante, Effects of Quenching and Partial Quenching on Penguin Matrix 
Elements, JHEP 0110 (2001) 037, [hep-lat/01 080 10]. 

[31] T. Bhattacharya et al, Calculating e'/ e using HYP staggered fermion, Nucl. Phys. B (Proc. Suppl.) 
140 (2005) 369-371, [hep-lat/0409046]. 

[32] T. Bhattacharya et al, Weak Matrix Elements for CP violation, Nucl. Phys. B (Proc. Suppl.) 106 
(2002) 311-313, [hep-lat/01 11004]. 

[33] RBC Collaboration, C. Aubin et al., Systematic effects of the quenched approximation on the strong 
penguin contribution to s'/e, Phys. Rev. D74 (2006) 034510, [hep-lat/0603025]. 

[34] M. Golterman and E. Pallante, Quenched penguins and the Al = 1/2 rule, Phys. Rev. D74 (2006) 
014509, [hep-lat/0602025]. 

[35] M. Golterman and S. Peris, Analytic estimates for penguin operators in quenched QCD, Phys. Rev. 
D68 (2003) 094506, [hep-lat/0306028]. 

[36] K. Choi and W. Lee, Re(A Q ), Re(A 2 ) and RG evolution for N F = 3, Nucl. Phys. B (Proc. Suppl.) 140 
(2005) 384-386, [hep-lat/0409020]. 

[37] L. Maiani and M. Testa, Final state interactions from Euclidean correlation functions., Phys. Lett. 
8245(1990)585. 

[38] M. Luscher and U. Wolff, How To Calculate The Elastic Scattering Matrix In Two-Dimensional 
Quantum Field Theories By Numerical Simulation, Nucl. Phys. B339 (1990) 222. 

[39] H.R. Fiebig et al., On the I = 2 channel % — % interaction in the chiral limit, Nucl. Phys. B (Proc. 
Suppl.) 73(1999) 252. 

[40] CP-PACS Collaboration, S. Aoki et al, I — 2 pion scattering phase shift with Wilson fermions, 
Phys. Rev. D 67 (2003) 014502, [hep-lat/0209124]. 

[41] Changhoan Kim, Al = 3/2 K — > n% with physical final state, Nucl. Phys. B (Proc. Suppl.) 140 (2005) 
381-383. 

[42] K. Rummukainen and S. Gottlieb, RESONANCE SCATTERING PHASE SHIFTS ON A NON-REST 
FRAME LATTICE, Nucl. Phys. B450 (1995) 397, [hep-lat/9503028]. 

[43] N. Christ et al., Finite Volume Corrections to the Two-Particle Decay of States with Non-Zero 
Momentum, Phys. Rev. D 72 (2005) 114506, [hep-lat/0507009]. 

[44] C. Kim et al., Finite-Volume Effects for Two-Hadron States in Moving Frames, Nucl. Phys. Will 
(2005) 218-243, [hep-lat/0507006]. 

[45] L. Lellouch and M. Luscher, Weak transition matrix elements from finite-volume correlation 
functions, Comm. Math. Phys. 219 (2001) 31-44, [hep-lat/0003023]. 

[46] RBC Collaboration, T. Yamazaki et al, Al = 3/2 K — > %% calculation with interaction effect on finite 
volume, PoS LAT2006 (2006) 100. 



28 



Progress in Kaon Physics on the Lattice 



Weonjong Lee 



[47] H. Leutwyler and M. Roos, Determination of the elements V us and V u( i of the Kobayashi-Maskawa 
matrix, Z. Phys. C 25 (1984) 91. 

[48] E. Blucher et ai, Status of the Cabibbo Angle, Proceedings of the CKM 2005 Workshop, 
[hep-ph/05 12039]. 

[49] J. Bijnens and P. Talavera, X/3 decays in chiral perturbation theory, Nucl. Phys. B669 (2003) 
341-362, [hep-ph/0303013]. 

[50] D. Becirevic et ai, The K —* n vector form factor at zero momentum transfer on the lattice, 
Nucl. Phys. B705 (2005) 339-362, [hep-ph/0403217]. 

[51] C. Dawson et ai, Vector form factor in K12 semileptonic decay with two flavors of dynamical 
domain-wall quarks, [hep-ph/0607162]. 

[52] JLQCD Collaboration, N. Tsutsui et ai, Kaon semileptonic decay form factors in two-flavor QCD, 
PoS LAT2005 (2005) 357, [hep-lat/0510068]. 

[53] M. Okamoto, B, D, K decays and CKM matrix from lattice QCD, FPCP04 proceeding, 
[hep-lat/04 12044]. 

[54] UKQCD Collaboration, PA. Boyle et ai, Lattice computation of the first moment of the Kaon 's 
distribution amplitude, [hep-lat/0607018]. 

[55] J.W. Chen and I.W. Stewart, Model Independent Results for SU(3) Violation in Light-Cone 
Distribution Functions, Phys. Rev. Lett. 92 (2004) 202001, [hep-ph/031 1285]. 

[56] QCDSF/UKQCD Collaboration, V.M. Braun et ai, Moments of pseudoscalar meson distribution 
amplitudes from the lattice, [hep-lat/0606012]. 



29 



